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Abstract 

A general first-principles theory of dissociative recombination is developed for highly-symmetric 
molecular ions and applied to H3O"*" and CHj]', which play an important role in astrophysical, 
combustion, and laboratory plasma environments. The theoretical cross-sections obtained for the 
dissociative recombination of the two ions are in good agreement with existing experimental data 
from storage ring experiments. 
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The dissociative recombination (DR) of H3O"'" and CH.^ ions with electrons is an impor- 
tant process in combustion, astrophysical, and laboratory plasma environment. For example, 
DR of HsO"^ is believed to be the major source for water production in interstellar clouds 
and tails of comets, while the CH3" ion triggers a series of reactions at the origin of formation 
of large hydrocarbons in the interstellar space. 

Description of DR of triatomic molecular ions with low-energy (<2 eV) electrons, de- 
veloped in recent years |2[, takes into account vibrational, rotational, and nuclear spin 
degrees of freedom. It provides detailed information about the recombination mechanism, 
about which rotational and vibrational states of the molecule participate, and about the 
final distribution of products. However, the full implementation of the method is very de- 
manding numerically, if not impossible for four-atomic molecules such as HsO^ or CH^, 
mainly due to the number of degrees of freedom involved. For Hjj' js^^ and HCO^ js 2|, a 
simpler approach has been developed, in which one retains only such ingredients that have 
been found in the detailed quantum-mechanical study to represent the main DR mechanism. 
As a result, the simpler approach has proven to be successful and to agree well with the 
detailed theory and experiment. 

The present study takes another step in the development of the low-energy DR theory by 
generalizing it to treat other highly-symmetric ions, and by applying it to HaO^ anc 
The method implemented is based on ideas developed in earlier DR studies BBS' 
it is significantly generalized and accounts at present for vibronic coupling that has been 
neglected in the earlier studies. This coupling is small in B.^ but appears to be important 
in ions with many electrons. As in Refs. 5|, l6[, the starting point of the approach is the set 
of ab initio potential surfaces of the neutral molecule and of the corresponding molecular 
ion, which are obtained numerically near the equilibrium position of the ion. 

Our study of DR in Hj]~ showed that at temperatures T >100 K the overall resonance- 
averaged DR rate is mainly determined by the capture of the incident electron into Rydberg 
states attached to excited vibrational levels of the molecular ion. Inclusion of the rotational 
structure at modestly higher temperatures does not appreciably change the overall DR 
rate. Therefore, in the simplified method we neglect rotation and assume that the excited 
vibrational modes are responsible for the electron capture. For molecules heavier than B.^, 
the approach should be valid even for temperatures less than 100 K. Another assumption in 
this approach is that predissociation is faster than autoionization once the electron has been 



CH+. 



a but 



2 



captured by the ion. This appears to be justified for the most of polyatomic ions 
With these assumptions, the DR cross-section a{Eei) for capture of the electron by an ion 
initially on its vibrational level v' via the temporary capture into a Rydberg state associated 
to the vibrational ionic level v is given by Eq. (14) of Ref. 

M^eO) = ^r..3, (1) 

where Egi = (^/c)^/(2me), k, e, and me are the kinetic energy, wave vector, charge, the mass 
of the incident electron, is the Bohr radius, r„ is the width of the autoionization resonance 
for the neutral molecule associated with the capture of the electron in the vibrational state 
V of the ionic core, u is the effective quantum number of the electronic wave function of the 
resonance. (Notice that symbols v and u are different. Also, an implicit assumption here 
is that a single electronic ionization threshold is relevant, whereby only one value of u is 
relevant.) Brackets () symbolically mean that the cross-section is averaged over the Rydberg 
series of resonances associated with the same vibrational state v of the ion jsl as well as over 
the rotational structure of the molecule. Below we omit the brackets for simplicity, but 
the displayed cross-sections below are assumed to be averaged over the Rydberg series. If 
there are several vibrational states contributing to the cross-section for the capture, then a 
corresponding sum over v should be evaluated. 

The above formula allows a quick estimation of the DR cross-section, if experimental 

fi n 

spectroscopic data (as in j^) or ab initio data (as in Refs. [2[r) exists that can be used 



to determine T^. As has been demonstrated in Refs. {2 5, 9|, this approach provides 
reasonable agreement with experiment for the molecular ions that have been considered to 
date. In order to calculate DR cross-sections for small and medium-size polyatomic ions 
with a modest computational ab initio effort, we rewrite the above formula in a different 
form. 

With the above assumptions about averaging over rotational Rydberg series, the DR 
cross-section a{Ef.i) for electron capture into vibrational level v and electronic state A, with 
the ion initially in vibrational level v' and the electron in state A', can be written in terms 
of the scattering matrix {v, A\S\v' , A') 

a=^^\{v,A\S\v',A')\\ (2) 

The matrix element in the above equation can be evaluated using the vibrational frame 
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{v,A\S\v',A') = j dQ{v\Q)SAMQ){QW)- (3) 
We assume that the ion is initially in the ground vibrational state \v') = |0). We use the 



with mode coordinates labeled 



normal mode approximation for the vibration (as in 
as Qi {Q = {qi,q2, ■ ■ ■}), the form S = exp(27ri/i) of the scattering matrix, and with an 
expansion of fiA,A'{Q) in a Taylor series around the equilibrium configuration Qq of the ion. 
This yields 

/iA,A'(Q) = /^A,A'(Qo) + V + . . . (4) 

Here we use dimensionless normal modes qi,q2, - ■■, which are related to the length- unit 
normal coordinates Si, S2, ■ ■ ■ as = Si jj^red^ / ^, where fired and uj are the reduced mass 
and the frequency of the normal mode. Retaining only zero- and first-order terms in the 
above expansion, Eq. ([2]) takes the form ^ 

-^iE.) = '4('-^Y\{vMm\'. (5) 



The index i in Vi and cTj is used to stress that the capture occurs into the qi vibrational 
mode excited by one vibrational quantum Vi. It is more convenient to use effective quantum 
numbers iy{Q) = n — fi{Q), where n is the principal quantum number, rather than quantum 
defects. In the harmonic oscillator approximation, the matrix element {vi\qi\0) = 6^-^i/\/2. 
This gives 

a,{E,i) = 4 (^!^Ye{hw, - E,i)gS.^,i . (6) 



k'^ \ dq, 

The fact that the matrix element {vi\qi\0) is nonzero only for Vi = 1 provides a rationalization 
for the general observation that electron capture typically occurs into the lowest vibrational 
level. Contributions ai{Ef.i) from all of the normal modes q^ must be added together to 
obtain the total cross-section. A particular contribution ai{Eei) should be set to (at 
this crudest level of approximation) once the energy E^i of the incident electron reaches 
the energy of one Vi quantum. This results in a sharp drop of the cross-section once the 
energy of the incident electron reaches the vibrational threshold - an effect that has been 
observed experimentally for several triatomic molecules. In Eq. this is accounted for by 
the Heaviside step function 9. We have also included the degeneracy factor g. 

The actual DR cross-section is, in general, smaller than the cross-section for electron 
capture obtained in the way described above. Autoionization following the capture should 



decrease it to some degree. Once the electron is captured, the branching ratio between pre- 
dissociation and autoionization probabihties determines the actual DR cross-section. The 
predisso elation probability is determined by the strength of coupling between the state |lj, A) 
formed in the collision with other states |i;-',A"). On the other hand, the "complex multi- 
channel resonance phenomenon" familiar from multichannel Rydberg spectroscopy (MQDT) 
0, [3 5 lead to an enhancement in the electron capture probability. In the previous 
detailed study of DR in B.^ it was found that the electron capture enhancement due to 
complex multichannel resonances compensates for the probability of autoionization. In par- 
ticular, it was found that the rotational-vibrational coupling helps to channel the energy of 
the captured electron towards higher-excited vibrational states of the ion. This leads to a 
decrease in the autoionization probability. In other words, even in B.^, the vibrational modes 
of the Hj]" -|- e~ system are coupled strongly enough to make the autoionization/dissociation 
branching ratio to be small. In larger polyatomic molecules, such intra-molecular vibrational 
energy transfer mediated by vibronic coupling should be even more efficient, although we 
have not made detailed calculations to confirm this in our applications to H3O''" and 0113". 
We anticipate that the upper-limit estimate of the DR cross-section developed in the present 
study should be close to the actual value. 

In order to use Eq. ([5]) we have determined the potential energy surfaces U{Q) of a 
few excited electronic states of the neutral molecule near the equilibrium position Qq of the 
ion. Then, the effective quantum numbers i^iQ) are obtained from the Rydberg formula 
U{Q) = U^{Q) — e^/(2aoz^^(Q)), U~^{Q) being the electronic energy of the ion. 

The energies U{Q) and U~^{Q) and the related quantum defects and effective quantum 
numbers, which are obtained from ab initio calculations, correspond to the adiabatic repre- 
sentation where fi and z> matrices are diagonal but depend strongly on Q, especially when Q 
is changing along one of the degenerate normal modes present in H3O"'" and CHj]": Because 
such modes couple electronic states of the e irreducible representation degenerate at Qo, 
the adiabatic approximation cannot be used in a straightforward manner. Fortunately, it is 
possible to construct a model Hamiltonian (similar to the one constructed in 14l4l6l|) for the 
electron-ion interaction in the diabatic representation. The same representation can then 
be used to construct the model matrices fi{Q) or O^Q) , which are not diagonal (as in any 
diabatic representation) and contain a few parameters that can be fitted to reproduce ab 
initio energies. 
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FIG. 1: H3O effective quantum numbers v{Q)- The two left panels show the dependence of u for 
several electronic states of H3O as a function of the x-component of the first doubly-degenerate 
vibrational mode 52 of HaO"^. The right panels show the same states for the second doubly- 
degenerate vibrational mode 54. For 7^ electronic states are characterized by Cg irreducible 
representations, i.e. a' or a" . 



The ah initio calculation has been made using the Columbus software suite [17|. In 
the calculation, the electronic wave function of the ionic core has been represented by the 
standard cc-pvtz 18|] basis. For an accurate description of the excited electronic states 
of H3O, s,p, and d diffuse universal Rydberg electronic orbitals optimized to represent 
the Rydberg wave functions have been added to the ionic core orbitals. The number 
of additional Rydberg basis functions and the level of the configuration interaction was 
increased until the quantum defects converged to the needed accuracy. The origin q2x = 
lix = corresponds to the equilibrium of HsO"*". 

The HaO"*" ion at equilibrium has a C-^v pyramidal form with three H atoms forming the 
equilateral base of the pyramid. It has four normal modes, and two of them are doubly 
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degenerate, namely the second q2 and the fourth ^4 modes (in order of increasing vibrational 
frequencies). We expect that the incident electron is captured by the ion due to the excitation 
of one of the degenerate modes, because the vibration along these modes breaks the C^y 
symmetry of the electronic states and, therefore, produces an additional vibronic coupling 
that is responsible for the large DR probability. Figure [T] shows the effective quantum 
numbers v{q2x) and ^{qix) on a limited interval of v for several electronic states obtained 
from H3O ah initio energies (calculations were performed up to z/ = 5). 
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TABLE I: The fitted coupling parameters. Lines of numerical values correspond to the triads of 
electronic states included in the present treatment. The e pairs of the triads are shown in Figs. 
[T] and O Frequencies uj are in cm~^ and calculated using the normal mode approximation and, 
therefore, somewhat different from the experimental values. In this study we enumerate vibrational 
modes in the order of increasing energy. Spectroscopic notations are different: The present (72 and 
(74 modes for both ions are modes 4 and 3 correspondingly. 

As one can see, three pairs of electronic states near v{Qq) = 3.7, 4.02, and 4.12 are 
degenerate at the C^^ geometry and are split by the Jahn- Teller coupling for distorted 
geometries. At Qo these electronic states corresponds to pir, dir and d6 states in the united 
atom approximation. Our initial intention was to obtain Jahn- Teller coupling parameters 
from the three pairs in order to obtain the DR cross-section as had been done in the previous 
study of Hj^. However, as is clear from the right lower panel of the figure, the Ade state is 
also coupled to 5sai. Therefore, we decided to include in the model Jahn- Teller Hamiltonian 
(and in the model matrix of effective quantum numbers), the corresponding coupling. The 
model matrix in the diabatic basis of two e and one ai electronic states can be written as 
a sum of unperturbed diagonal matrix i>'-°\ which is q2- and g4-independent, and the matrix 
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FIG. 2: Comparison of the present theoretica. 
plotted experimental data are from Refs. 



and experimental H3O''' DR cross-sections. The 



2[h22|. The theoretical curves include a convolution 



with an experimental anisotropic distribution E\\ =0.1 meV and E_\_ = 2 meV of electron energies. 
The theoretical curves also accounts for the toroidal correction as discussed in [23], which has a 



minor effect on the cross-section in the absence of autoionizing resonances. The arrows indicate 
threshold energies for the q2 and 54 modes. 

of perturbation h'^^\ which is linear in q2 and q^. The perturbation is caused by the 
symmetry breaking. Diagonal elements of u^^^ are zero. Hence the total matrix u can be 
written as 



,(0) 



-up 



(0) 



f^eaP^ 
,,(0) 



lip 



(7) 



where p and if are polar versions of the x- and y- components of q2 or q^. The parameters 
Kee and Kea are real, p- and independent, and obtained by a fit to the ab initio results 
shown in Fig. [TJ The eigenvalues of the above matrix should be equal to the values obtained 
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FIG. 3: Effective quantum numbers and DR cross-sections (theoretical and experimental [20[) for 
CHg". Figure shows only a limited interval of v{qx) values to demonstrate that in the CHj^ case, 
the coupling between e' and a'^ (D^h) electronic states is strong at distorted geometry because the 
quantum defects of the e' and a'^ electronic states at Qo a-re almost degenerate. For ^ the 
electronic states are characterized by the C2v irreducible representations, i.e. oi, 02, 61, or 62. 

from ab initio energies. The fitted parameters K^e and n^a are given in Table [I] for both ions. 

Using the fitted values Kee and Kga, the contribution to the total DR cross-section from 
each triad of the e and Oi states for each vibrational mode is given by Eq. The 
degeneracy factor has the value g = 2 ioi the e ~ e interaction: Only two projections (say 
the piT"^, dn^, or d6^ electronic states in the united-atom approximation) of the angular 
momentum of the incident electron on the molecular symmetry axis can excite vibration 
of the ionic core. For the e ~ a interaction g = 4: If the incident electron arrives in an 
e electronic state included in the matrix of Eq. ([7]), it can excite only one component of 
the doubly degenerate q2 or modes by the e ~ a coupling. There are two such electronic 
states e. However, if the electron arrives in the a state included in the matrix of Eq. ([7]), 
both vibrational components of the q2 or q^ modes can be excited. The a state is coupled 
to the two e electronic states and to both vibrational components of q2 or ^4 with the value 
g = A. 

Figure [2] shows the resulting DR cross-section for H3O''" and compares it with data from 
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storage ring experiments. To compare with the 
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experimental data, the non-uniform experimental distribution of collision energies between 
ions and electrons has been averaged over for the theoretical result. We have also made a 
similar calculation for the CB.f ion, which has D^h symmetry at the equilibrium and can be 
treated in the same way as HsO^: It has two E' doubly degenerate modes that contribute 
to the DR cross-section. The result of the calculation is summarized in Fig. [31 Here 
the theoretical cross-section also agrees well with the experimental data except the energy 
interval (0.2-0.4 eV) between the two vibrational thresholds. The disagreement could be 
explained by a presumably high rotational temperature in the experiment. The rotational 
structure is not accounted for in the theoretical curve. 

In conclusion, we stress that the present approach for theoretical determination of the 
DR cross-section in small polyatomic closed-shell ions at energies below 1 eV is general, 
simple in application, and gives good agreement with experiment as well as with the detailed 
theoretical study made for the benchmark problem of DR in H^. 
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CH.^ and T. Oka for helpful comments on the manuscript. This work has been supported 
by the National Science Foundation and in part by the Department of Energy, Office of 
Science. VK also acknowledges support from the RTRA network Triangle de la Physique. 
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